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Abstract 

We propose an extended quantum mechanical formalism that is based on a 
wave operator which is related to the ordinary density matrix via p = qqK 
This formalism allows a (generalized) unitary evolution between pure and 
mixed states. It also preserves much of the connection between symmetries 
and conservation laws. The new formalism is illustrated for the case of a 
two level system. 
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Several proposals motivated by various considerations for generalizing the quantum 
mechanical formalism have been made todate. In these programs one disposes of a fun- 
damental quantum mechanical principle such as linearity, locality, or unitarity. Weinberg 
suggested a nonlinear generalization and proposed precision tests of nonlinear corrections 
to quantum mechanics [|l|. Motivated by the apparent breakdown of unitarity in the black- 
hole evaporation process. Hawking proposed that a synthesis of quantum mechanics and 
general relativity requires giving up unitarity , and to some extent locality 0] . A model 
which gives up both properties, was constructed by Marinov 0]. As a linear and local 
phenomenological implementation of Hawking's proposal, Ellis, Hagelin, Nanopoulos and 
Srednicki (EHNS) |^, and Banks, Peskin and Susskind (BPS) 0], suggested a modified 
Liouville equation for the density matrix p. 

In particular, BPS showed that the requirements of linearity, locality in time, and 
conservation of probabilities, lead to a modified equation with a "generic form" : 

ihdtp = [H,p]+iY^ KmiQmQnP + pQmQn " 2Q„,pQm), (1) 
71, m 

Here, Qn are any Hermitian operators and hnm is c-number hermitian matrix. A sufficient 
but not necessary condition ensuring the positivity of p, is that the matrix h is positive. 
Equation (|l|) does not preserve trp^. Thus pure states can indeed evolve to mixed states 
i- 

Similar equations can be obtained from ordinary quantum mechanics for a subsystem 
interacting with an environment Nevertheless, when gravity is involved one can argue 
that the relevant "micro-environment" is hidden by black-hole horizons and is in principle 
unobservable. This would render equation (|T|) a fundamental modification of quantum 
mechanics, rather then an artifact of interacting with an environment. 

Modified evolutions like ([^) were applied in various cases. EHNS proposed that the cor- 
rections induced might be observed in the ultra-sensitive Kq — Kq system. Furthermore, 
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Ellis et. al. |^ and Huet and Peskin [|T^] examined the possibility that the observed CP 
violation in the Kq — Kq system is, to some extent, due to non-quantum mechanical cor- 
rections. Related modifications where also proposed in connection with the "measurement 



problem" , in order to generate a von- Neumann reduction for macroscopic systems [0, |T2 

In what follows, we propose a different approach. It is also based on the Liouville 
equation but not for the ordinary density matrix. It constitutes a linear, local, and unitary 
extension of quantum mechanics. To this end, consider the density matrix in ordinary 
quantum mechanics and focus first on the case of a pure state. By analogy with the 
relation p = \ip){ip\, let us define the operator g by |13 



P = Q6^- (2) 

If g satisfies a Liouville equation 

ihdtg = [H, g], (3) 

it is easy to see that the density matrix p = gg^ also satisfies a Liouville equation with 
the same Hamiltonian. The initial condition may be specified in terms of the "square root 
operator" g, rather than p. Thus, if the system is determined at t = to by an ordinary 
complete set of measurements to be in the state {ipo), (or p = \4'o){4'o\), this sets the initial 
condition for equation (|): 

^(t = to) = |^o)(^o|. (4) 

Now we observe that and imply p(t = to) = g{t = to), and since both quantities obey 
the same equation of motion this relation holds at any subsequent time. The expectation 
values of any observable A is obtained by the standard expression: 

tr Ap trAg 

= - — = — (5) 
trp tr^ 



Hence eqs. (2-5) are equivalent to ordinary quantum mechanics [14 



It is therefore interesting to question whether eq. (3) can now be used as a new starting 
point for a quantum mechanical extension. We shall assume that q is from now on a 
general operator (not necessarily a projector) still obeying the initial condition (4), and 
that expectation values are still obtained by the standard expression 



The density matrix however is from now on obtained via p = qq\ 

The hermiticity and positivity of p is automatically ensured by p = qq^ . We need that 
the modified equation conserve probabilities, i.e, dtiip = dttigg'^ = 0, but not necessarily 
purity. The most general linear and local generalization of eq. (|^) which satisfies this 
condition can be written as: 



Here, L, R, Ki and Kj are any Hermitian operators, Qij are real coefficients, and the 
summation convention was used. 

Eq. implies that the density matrix obeys: 



The "primary" object g can not be eliminated from eq. (|^) which therefore cannot be 
rephrased in terms of p only. Thus unlike the case of eq. (1), p plays here the role of 
a "secondary" object. Eq. (§) also indicates that the term Lg in eq. gives rise to a 
redefinition of the Hamiltonian and that the term gR can be eliminated. Indeed, the gauge 
transformation, g — > gU, where U is a unitary operator, does not affect expectation values 
and can be used to recast eq. (j^) into the form: 



ihdtg = [H, g] + Lg + gR + g^KigR'^. 



(7) 



ihdtp = [H + L,p]+ gij{KigKjg^ - h.c). 




ihdtg = Hg + gijKigK. 



(9) 



3 



-ij\R-H)dt' 



. Without the last term 



where, H = H + L, Kj = UKp-\ and U = exp 
this is simply a Schrodinger-like equation for the operator g. 

To further analyze eq. (0) we construct a Hilbert space. It is defined as the linear space 
C = {g} of solutions of eq. with all possible initial conditions at any to. With the inner 
product defined as: 

(^1, ^2) = tr^l^2, (10) 

C becomes a Hilbert space. It follows from eq. that this inner product is conserved and 
hence the generalized dynamics suggested here manifests in £ as a unitary evolution. The 
inner product (|1^) may be regarded as an extension of the ordinary quantum mechanical 
inner product. If the corrections induced after t = to by the new terms in the evolution 
eq. are small, (^1,^2) — KV^i|V'2)P- Note also that expression (j^) for the expectation 
value of of an observable A can be now re-expressed as: 

(A) = (11) 
Equations (|^J|) , and ( p!0| - |TT|) suggest that g should be interpreted as a generalized "wave 



operator". The new feature here however, is that trp^ = tr(^^^)^ is not conserved. This 
manifests the new aspects of our unitary evolution as transition between pure and mixed 
density matrices (p). 

The generalized unitarity, namely the conservation of the inner product (0), can be 
clarified by rewriting eq. (|^ in the Hilbert space C. For simplicity let us consider a system 
with a finite, A^-dimensional, Hilbert space and perform the extension described above. 
The extended, A^^ dimensional, Hilbert space C can be spanned by a hermitian basis of 
A^^ — 1 SU(N) matrices and the unit operator: 

Q = -^{Qo'I + QiTi), (12) 
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where Tj are SU(N) generators and Qa are A^^ complex numbers. In this basis, the gener- 
ahzed inner product between any two solutions is given by an ordinary vector product in 
a A^^-dimensional Hilbert space: 

(^1, ^2) = QlaQ2a- (13) 

a=0 

We can also express eq. (J^) in this basis as a Schrodinger-like equation: 

iUdtQa = KabQb = (^Ir^ + SHab)Qb- (14) 

The condition for conservation of probabilities (and unitarity) is simply that the gener- 
alized Hamiltonian, Tiab, is hermitian. The deceptive similarity of eq. (|I4|) and ordinary 
quantum mechanics Schrodinger equation in an N"^ dimensional space notwithstanding, we 
emphasis that the only relevant, physical degrees of freedom are in those of the original 
(A^ dimensional) Hilbert space. 

Next we would like to express the observables Ai as hermitian operators in C. In general 
we have in C independent hermitian operators. Therefore the mapping 

Ai^Ae Oc (15) 

of the original (A^^) observables Ai into the set of hermitian operators Oc in £^ is not one 
to one. This mapping is constrained by demanding that 

tlAiQQ^ = H H Qa{A)abQb, (16) 
a=0 6=0 

i.e., that {A) is expressible in £ as a "standard" expectation value with respect to the 
"amplitudes" Qa- We also require that the mapping (|T3|) preserves commutation relations. 
Therefore, an A^-dimensional representation of SU(N) is mapped into an A^^ dimensional 
representation of SU(N) in Oc-, Ti —>■ %. The linear transformation maps a general ob- 
servable Ai = Cjol + CiaTa to Ai = Ci^X + J^a^iaXi- The Operator Ai G Oc still has the 
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same eigenvalues as the original operator Ai. However all the eigenvalues are now A^-fold 
degenerate. Another set of operators, denoted by Vj, which remove the degeneracy of Ai 
do not correspond to observables. It can be shown that the role of the new terms in eq. 
(0) or 6H in eq. ( p!4D is to generate correlations between Ai and Vj, which in turn induces 
the transition to a mixed density matrix. 

It was noted by Gross ||16[ and by Ellis et. al. [^j, that linear modifications of the 
evolution laws for the density matrix (e.g. eq. ([l|)) generally breaks the one to one corre- 
spondence between symmetries and conservation laws. We now show that in the present 
formalism, this correspondence is partially restored. An observable A G Oc that is a 
constant of motion satisfies: [.A, 7Y] = 0. Hence the unitary operator T = exp{—ieA/h) 
commutes with the unitary evolution operator U = exp{—itT-l/h), and A generates a sym- 
metry in C. The converse is not generally true. Since C is A^^-dimensional, not all the 
hermitian operators in Oc may be mapped back to hermitian operators in the original A^- 
dimensional Hilbert space. Therefore, if some hermitian operator Q generates a symmetry 
in C and its expectation value, glit)QabQbit), is conserved, it still may not correspond to 
an observable. 

To illustrate the general discussion above let us consider as an example the simple two 
level system (e.g. a spin half particle in a constant magnetic field). The mapping between 
the original 2-d Hilbert space and the 4-d Hilbert space C will be spelled out in detail. Let 
the "free" Hamiltonian be given by 



H = Eq -{ — huja^. 



(17) 



We have seen that the terms Lq and qR in eq. (|^ can be absorbed by a redefinition of H 
and K'-. Therefore, the modified eq. will be taken as: 



ihdtQ = [H, g] + KqK', 



(18) 



6 



where K and K' are functions of the Pauh matrices, and will be assumed to be time 
independent. Energy conservation, -§^{H) = ^ ^ (ggf = 0, implies that [crs,-^'] = 0, hence 
K = a-s- This leaves three unknown parameters which determine K': 



K' = aoi + 13 02 + Acts. 



(19) 



When reexpressed in the four dimensional Hilbert space C the modified dynamics cor- 
responds to eq. with: 

/A il3 -ia 0\ 

—i/3 —A a 

+ia -X P 

V a f3 XJ 
The observables in this model are combinations of cTj and the unit operator. The 

mapping, (jj — > 5j G Oc is: 



6n 



(20) 



1 

— (7 k 

2 



^abk ) • 



(21) 



The Si are a 4-dimensional representation of SU{2), preserving the commutation relation 



= ieijkSk- The mapping (|2l| ) was constructed so as to satisfy eq. ([T6|). The 
operators, "Dj, which remove the degeneracy of Si have also been explicitly constructed. 
The latter indeed do not correspond to observables. 

It can now be verified that 1S3 is a constant of motion, i.e. [1S3, TY''"^ + dTi] = 0. We also 
notice that since the energy operator, Eol + huS^, is not the mapped original Hamiltonian: 



ab 



ihueabs, '^'"^ does not correspond to an observable in C. 



The present model differs qualitatively from the model of BPS or EHNS: while eq. 
(|I|) yields in general exponentially decaying (or exponentially increasing) solutions, our 
modifications are oscillatory. Indeed the general solution of equation (|T^ is 



6a= J2 c^,gi,aexp{-iXat), 

fj.=0,3 



(22) 



where Aq, and Qaa are the (real) eigenvalues and eigenvectors, respectively, of Hab- 

As an example consider the special case where only A in eq. (|2y) is non-vanishing, and 
the spin is found at t = to iii the state \ipo) = cos || If) + sin || l^). The solution in this 
case is given by 

m = { , ■ ; ^ : ... i- (23) 



iXt 



i sin(r7)e^(^+^)* sin^(r//2)e- 
The resulting density matrix, p = gg^ , oscillates periodically between a pure and mixed 



state. For example in the simple case rj = 7t/2 



-iu)t 



Pit) 



and trp2 = ^ + ^ cos2(2At). 



cos(2At) ■ 



2 Ve*^*cos(2At) 



(24) 



Observable effects due to these modifications can in principle be searched for in neutron 



interferometry experiments |T^. In such interference experiments, one typically measures 
an observable of the form 

m = h 1, (25) 



1 e'" 

where 9 is determined by the experimental set up. The expectation value of A is given in 
our case by 



{A) 



1 + sinr] cos^(r7/2) cos((c<; + 2A)t + 9)+ sin^{ri/2) cos((c<; - 2A)t + 9) 



(26) 



The correction is indeed oscillatory. This should be contrasted with the exponential 
exp(— 2AsAri/5t) decay of the interference obtained by EHNS. 

What are the present experimental bounds pertinent to the three new parameters of the 
two level system? We can use the two slit experiments of Zeilinger et. al. [|18| with a 20^4° 



neutron beam, and the analysis of Pearle |jT9|, to constrains the generic parameter A to 
A ~ 10 (sec) ~ 10^^'^ GeV. The constrains of the same experiment on the corresponding 
parameters in the EHNS model is ~ 100 times stronger (~ 10^^^ GeV). The exponential 
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factor modifies the interference contrast during the short flight time (to — 10^^ sec) by 
(1 — 2\EHNsto)- In the present case the extra oscillation can be subsumed into slow 
"beating" ~ cos(2Ato) — (1 — 2(Ato)^), causing a much weaker reduction of the contrast in 
the interference pattern. 

We found that our modification induces Kl Kg mixing generating CP violation in the 
two level Kq — Kq system in a similar fashion as in the EHNS model. However, this mixing 
predicts a phase of the CP violating parameter e, which is of by 7r/2 just as in the case 



of the EHNS model |]10l- Hence our modification can account only for a small part of the 
CP violation observed in the Kq — Kq system. This leads to the generic upper bound of 
order ~ 10^^^ Gev, of the same order as Mj^/Mpi which could be expected on dimensional 
grounds if CP and/or CPT violations are due to effects of quantum gravity. The hundred 
fold larger parameter allowed by the neutron interference experiment in our model could be 
important. In particular this renders smaller yet experimentally detectable CPT violations, 
more likely in the present framework. 

We have constructed a formalism based on an operator generalization of the wave 
function which is linear, local, and unitary. As a consistency check of this proposal we note 
that to some extent the proposed formalism can be embedded in the framework of ordinary 
quantum mechanics. We can interpret qq^ and q as the reduced density matrices of a sub- 
system and an environment, respectively. The generalized Hamiltonian Tiab in eq. (|T4|) and 
the amplitudes Qa can then be interpreted as the Hamiltonian and wave function of the total 
system, while the new terms in eq. or (0) as describing an interaction between the 

sub-system and the environment. Therefore, the consistency of the proposed equation of 
motion follows from quantum mechanics. Nevertheless, postulate (4), git = to) = \'^o){'^o\, 
which sets the initial condition for eq. goes beyond any ordinary quantum mechanical 
scheme. It would amount in quantum mechanics to an additional requirement that after 
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carrying a complete set of measurements on the sub-system, the wave function of the 
environment becomes identical to that of the system. This additional constraint is not 
satisfied in quantum mechanics. Therefore the predictions of this formalism will generally 



differ from that of a quantum mechanical system with an environment |20 



Finally, we note that the proposed formalism may also be relevant to the information 
problem in black hole evaporation and to the measurement problem. In the latter case, for 
large systems the modified evolution might under appropriate conditions give rise to loss 
of coherence which amounts to a measurement. 
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